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Abstract. We present the results of a nonadiabatic, linear stability analysis of models of very low-mass stars 
(VLMSs) and brown dwarfs (BDs) during the deuterium burning phase in the center. We find unstable fundamental 
modes with periods varying between ~5 hr for a 0.1 Mq star and ~1 hr for a 0.02 Mq BD. The growth time of 
the instability decreases with decreasing mass and remains well below the deuterium burning time scale in the 
mass range considered (0.1-0.02 Mq). These results are robust against variations of the relevant input physics 
in the evolutionary models. We identify possible candidates for pulsational variability among known VLMSs and 
BDs in nearby star forming regions whose location in the HR diagram falls within or close to the boundary of the 
instability strip. Finally, we discuss the possibility that the variability observed in a few objects with periods of 
~1 hr can be interpreted in terms of pulsation. 
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1. Introduction 

Stellar pulsation represents a unique probe to study the 
internal structure and evolution of stars. This technique 
has been applied to many classes of stars, both on the 
main sequence and in earlier/later stages of their evolu- 
tion (e.g., Gautschy & Saio 1995, Marconi & Palla 1998). 
New classes of pulsators are being discovered across the 
HR diagram from ground-based observations and with the 
advent of space missions (MONS, COROT) the potential 
of discovering solar-like oscillations in other stars promises 
a major improvement in our understanding of stellar in- 
teriors (e.g., Thompson, Cunha & Monteiro 2003). 

In this Letter, we propose the identification of a new 
class of pulsators, VLMSs and BDs, in which the in- 
stability is induced by deuterium (D-) burning in their 
centers. Due to the high sensitivity of nuclear reactions 
on temperature, fully convective stars can become pul- 
sationally unstable during the burning phases. Gabriel 
(1964) was the first to suggest that such mechanism, called 
e-mechanism, would make low-mass stars on the main 
sequence vibrationally unstable. Later on, Toma (1972) 
showed that pre-main-sequence (PMS) stars in the range 
0.2-2.0 Mq could also become pulsationally unstable dur- 
ing the D-burning phase and suggested a possible rela- 
tion of this phenomenon with the observed variability of 
T Tauri stars. In both cases, however, and for different rea- 
sons, the suggestion has not met with success and to date 
there is no observational evidence for the e-mechanism 



in any class of pulsating star where the destabilization 
is mainly induced by opacity effects in ionization regions 
near the stellar surface (the K-mechanism) . 

Indeed, no periodic variability, which could be inter- 
preted as pulsation, has been reported up to now in main 
sequence low-mass stars, which may not be surprising, 
given the weak dependence of the pp-chains on temper- 
ature and the resulting long e-folding times for the e- 
mechanism. For young PMS stars, the main difficulty to 
confirm observationally the results of Toma (1972) is due 
to the very short D-burning time scale, < 1 Myr, for 
stars with mass > 0.2 Mq. Moreover, the results obtained 
by Toma (1972) are optimised, in terms of D-burning 
time scale and e-folding time for the pulsation, since he 
adopted an initial D-abundance about ten times higher 
than the currently measured interstellar value and ignored 
D-deplction during protostellar accretion that can sub- 
stantially reduce its abundance at the beginning of the 
PMS phase of low-mass stars. As we will show below, 
VLMSs (<0.1 M Q ) and BDs in the earliest evolutionary 
stages have the appropriate physical conditions to circum- 
vent the limitations of previous attempts. 

2. Instability induced by D-burning in VLMSs and 
BDs 

BDs are expected to begin the PMS phase with convective 
interiors and with the full amount of interstellar deuterium 
since during the preceding phase of protostellar accretion 
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their centers are too cold to start any nuclear reaction. 
The ignition of deuterium requires temperatures around 
10 6 K that can only be reached in protostars more massive 
than ~0. 1-0.2 Mq, depending on the mass accretion rate 
(Stahler 1988). Thus, very low-mass objects need to con- 
tract somewhat in the PMS phase before deuterium can 
start burning and being depleted. Once the critical tem- 
perature is achieved, the D-burning phase occurs on a time 
scale that varies between td ~2.5 Myr for a 0.1 Mq star 
and td ^20 Myr for a 0.02 Mq brown dwarf (Chabrier et 
al. 2000). 

Since the energy generation rate, c nuC7 scales approxi- 
mately with the 12-th power of the temperature, the rate 
of combustion is slower for the lowest mass objects. The 
same temperature sensitivity of e nuc is at the root of the 
instability induced by any temperature variation: since 
Se nuc /e nuc ~ 12 ST/T, a small T-perturbation induces a 
variation of e nuc which is an order of magnitude bigger. 
(Note that e nuc is also linearly proportional to density, but 
no destabilizing effect is expected from variations in this 
quantity). In terms of pulsation analysis, the time scale for 
the growth of the instability, T growt h, is inversely propor- 
tional to 5e nuc and should be shorter than the D-burning 
time for the e-mechanism to operate. Therefore, in order 
to test its viability, it is important to follow numerically 
the growth of the instability in time since the onset of the 
D-burning phase. 

3. Evolutionary models and pulsation analysis 

Our analysis focuses on objects in the mass range 0.1-0.02 
Mq, close to the minimum value for D-burning (0.013 Mq, 
Chabrier et al. 2000). The evolutionary models are based 
on the input physics described in Baraffe et al. (1998) and 
on the so-called "NextGen" grainless atmosphere mod- 
els (Hauschildt et al. 1999). Models are calculated for a 
solar metallicity, with an initial deuterium mass fraction 
Do = 2 x 10~ 5 , characteristic of the local interstellar 
medium (Linsky 1998), and with mixing length parame- 
ter equal to the pressure scale height. During the early 
phase of D-burning, brown dwarfs are relatively hot, with 
T e ff > 2300X, and their atmospheric properties are still 
well described by dust-free models. 

3.1. Linear stability analysis results 

We performed a non-adiabatic linear stability analysis 
of the model structures during the D-burning phase, 
searching for the presence of unstable radial eigen modes. 
Unstable modes characterize oscillations around the hy- 
drostatic equilibrium configuration. If perturbations have 
time to grow, they could reach large amplitudes and re- 
sult in periodic phases of expansion and contraction, with 
a pulsation period related to the dynamical time scale of 
the object Tdyn ~ (Gp) -1 / 2 . References on the pulsation 
code used for the present purpose can be found in Lee 
(1985) and Baraffe et al. (2001). In all models analyzed, we 
find unstable fundamental modes, the results being rather 



robust against variations of the relevant input physics (see 
Sect. 3.2). We also find unstable first overtones, but their 
stability properties are sensitive to the input physics and 
their growth rate is usually smaller than that of the funda- 
mental mode. Higher overtones remain stable during the 
whole D-burning phase. We thus restrict the rest of our 
study to fundamental modes. 

In the mass range of interest, fundamental mode peri- 
ods Pq vary between ~lhr and ~5 hr, close to Td yn - For 
a given mass, Pq decreases during the D-burning phase, 
since contraction proceeds and p increases. This is illus- 
trated in Table 1 which summarizes the characteristics 
of VLMSs and BDs and their main pulsation properties 
during the D-burning phase for different masses. We de- 
fine Tg r owth as the e-folding time characterizing the growth 
time scale of the pulsation amplitude for the fundamental 
mode. Fig.^displays Pq and T growt h for a 0.1 Mq VLMS, 
and 0.06 and 0.03 Mq BDs. Figure Hand Table 1 show 
that at the beginning of the D-burning phase, r grow th is 
shorter, by at least a factor 10, than the D-burning time 
scale r D , which is 2.5, 4.2 and 9.3 Myr for 0.1, 0.06 and 
0.03 M Q , respectively. Note that only unstable models are 
shown in Fig. since the structure becomes stable again 
before complete depletion of D. It remains small during 
a significant fraction of the D-burning phase in the low- 
est mass BDs, indicating that the perturbations may have 
time to grow to reach significant amplitudes. The loca- 
tion of the D-instability strip in the HR diagram is shown 
in the top panel of Fig. along with evolutionary tracks 
of selected masses and isochrones. The curves of constant 
period cut the strip almost horizontally (nearly constant 
luminosity) and the shortest periods are found for the least 
massive BDs. Note that our models can also provide the 
locus of the D-instability strip in the color-magnitude di- 
agram (CMD) for different bands for comparison with ob- 
servations. We plan to present such diagrams in a forth- 
coming paper. 

Inspection of the differential work integral dW/dm 
throughout the structure of D-burning objects reveals the 
exciting (positive dW/dm) and the damping (negative 
dW/dm) zones. This quantity is plotted in Fig. for the 
0.03 Mq BD at the beginning of the D-burning phase. The 
positive work in the central region is due to the pertur- 
bation of the nuclear energy generation rate and charac- 
terizes the e-mechanism. As a rough approximation, the 
work integral in the central region follows 

W cx e nuc e T , (1) 

with e nuc the nuclear energy generation rate and 6t = 
( * n £ ™ ) p . The source of nuclear energy is due to the 
2 H(p,j) 3 He reaction with e^ ~ 12 — 14 in our objects, 
at temperatures T « 10 6 K. It is the combination of 
such high sensitivity of the nuclear energy generation rate 
with temperature, and the fully convective structure of 
low-mass objects, which yields a less centrally condensed 
structure, that allows the excitation by the e-mechanism. 
Indeed, VLMSs and BDs are characterized by a low cen- 
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Fig. 1. Evolution of the D-abundance (in units of the 
initial abundance Do = 2 x 10~ 5 in mass fraction), the 
fundamental mode period Pq (in hr) and growth time scale 
Tgrowth (in Myr) as a function of time during the D-burning 
phase of a 0.1 M© star (dashed lines), and a 0.06 M© 
(dash-dotted line) and 0.03 M© (solid line) BD. Note that 
only unstable models are shown. 



tral condensation with a ratio of central to mean density 
p c /p w 4.5 — 6. This ratio decreases with mass and for 
a given mass it remains nearly constant with age. The e- 
mechanism is usually unimportant for radial pulsation in 
stars because of the small amplitudes that can develop 
in dense central regions where nuclear burning occurs. 
However, stars less centrally condensed are thought to be 
potential sites for the e- mechanism (Gabriel 1964). Finite 
amplitudes of the eigenfunctions can indeed be expected 
in the central regions, as illustrated by the behavior of the 
pulsation amplitude Sr/r displayed in Fig.|3 

3.2. Uncertainties of the results 

In order to analyze the sensitivity of the present in- 
stability to input physics, we have computed evolution- 
ary models using atmosphere models based on different 
molecular linelists for TiO and H2O (see Chabrier et al. 
2000 for more details). The main effect is a change of 
the effective temperature (by ^100 K), and of the D- 
burning time scale. As a test on the initial conditions, 
we have constructed evolutionary sequences starting at 
temperatures T c < 5 x 10 5 K before central D-ignition. 
Moreover, we have explored the effect of the mixing length 
(imix = 1 — 2 Hp), since it can affect the early evolution 
of VLMSs and BDs (see Baraffe et al. 2002). All these 
variations yield only small quantitative differences, mainly 
in the D-burning lifetime and the pulsation growth time 
scale. However, the instability of the fundamental mode 
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Fig. 2. Top panel: Location of the D-instability strip 
(shaded area) in the HR diagram. Tracks for different 
masses and isochrones are indicated as labeled. The isope- 
riod curves (in hr) are shown by the thick lines within the 
strip. Bottom panel: Distribution of known VLMSs and 
BDs in Taurus (circles) and Cha I (triangles). Data points 
for Taurus are from Briceho et al. (2002), Luhman et al. 
(2003), and Luhman (2004), and for Cha I from Comeron 
et al. (2000). 

survives in all cases, with r grow th <C td during most of the 
D-burning phase. Finally, adopting for the 2 H(p,j) 3 He 
reaction the rate of Caughlan & Fowler (1988) or that 
of the NACRE compilation (Angulo et al. 1999) has no 
consequence. 

An important uncertainty stems from the treatment of 
convection, which is assumed frozen in the linear stability 
analysis, having neglected the perturbation of the convec- 
tive flux in the linearized energy equation. Such standard 
approximation is only justified by the argument that the 
convective time scale in the inner structure of our objects 
(Tconv ~ 10 6 — 10 s s ) is significantly larger than the pulsa- 
tion periods ( <J 10 4 s). Analysis of the effect of convection 
on the perturbation is beyond the scope of this study and 
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Fig. 3. Differential work integral dW/dm, in arbitrary 
units (solid line), as a function of interior mass (in units 
of the total mass) for a 0.03 Mq BD at the beginning of 
the D-burning phase (D-abundance D/Dq=0.71, age=2.4 
Myr). The pulsation amplitude Sr/r (dotted line), arbi- 
trarily normalized to 1 at the surface, and the nuclear 
energy generation rate e nuc (in erg/g/s, dashed line) are 
also displayed. 

in any case can only be qualitative given our poor knowl- 
edge of the interaction between convection and pulsation. 

4. Discussion 

Our results show that the perturbations excited in the 
center of D-burning VLMSs and BDs have smaller growth 
time scale than that of the D-burning phase and may thus 
have time to grow to the level of observable amplitudes. 
Unfortunately, our linear models cannot be used to es- 
timate the amplitude of the pulsations. More theoretical 
work is required and we hope that these initial results 
will stimulate the development of non-linear calculations 
in order to estimate pulsation amplitudes and of fully hy- 
drodynamical models to construct synthetic light curves 
for comparison with observations. In any case, it is impor- 
tant to realize that it may be difficult to distinguish pul- 
sations from other sources of photometric variability, such 
as rotation and spots or atmospheric events, that occur in 
VLMSs and BDs on similar time scales. Fortunately, nu- 
merous candidates exist which could provide good targets 
for detailed observational studies. 

As an illustration, the lower panel of Fig. [21 displays 
the location in the HR diagram of 25 spectroscopically 
confirmed VLMSs and BDs drawn from larger surveys in 
Taurus and Chamaclcon I (Briceho et al. 2002, Comeron 
et al. 2000, Luhman et al. 2003, Luhman 2004). The rather 
large error bars indicate the current uncertainties in the 
conversion from spectral types to effective temperatures 



(see discussion in, e.g., Luhman 2004), but the overall dis- 
tribution closely matches the predicted position of the in- 
stability strip down to the lowest masses with the shortest 
pulsational periods. 

Evidence for photometric variability of young VLMSs 
and BDs both in the field and in young clusters has 
been already obtained by several groups in the last few 
years (Bailer- Jones & Mundt 2001, Joergens et al. 2003, 
Zapatero Osorio et al. 2003). The amplitude of the ob- 
served variations vary between tens of mmag in the optical 
to —0.05-0.2 mag in the near-IR where these cool objects 
emit most of their energy. In several cases, periodic vari- 
ability has been reported with tentative periods in the 
range from half an hour to several hours. Very interesting 
results have been found in the a Orionis clusters (age —3 
Myr) where two objects (S Ori 27 and S Ori 28) have pe- 
riods around 3 hr and even shorter in the case of S Ori 31 
(-2 hr) and S Ori 45 (-0.75 hr) (Bailer-Jones & Mundt 
2001, Zapatero-Osorio et al. 2003). The minimum funda- 
mental periods that we find are close to ~ lhr, correspond- 
ing to the lowest mass studied (cf. Table 1). Although the 
period of S Ori 45 is slightly shorter than 1 hr, it is com- 
pelling that this object has precisely a mass which could 
be as low as 0.02 M Q (Bejar et al. 1999). Several other 
candidates in the e Orionis cluster appear to share similar 
properties (Scholz & Eisloffel 2004). So far, the observed 
variability has been interpreted in terms of rotation peri- 
ods, presence of photospheric cool or hot spots, interac- 
tion with accretion disks, and/or atmospheric events. For 
objects with the shortest periods (—1 hr), the inferred ro- 
tational velocities would exceed —100 km s _1 at or above 
breakup (e.g., Joergens et al. 2003). The pulsational in- 
stability induced by D-burning thus offers a new interpre- 
tation that relics on fundamental (sub)stellar properties. 
It also has the potential of providing direct information 
on the otherwise inaccesible internal structure of VLMSs 
and BDs. More fundamentally, it could provide the first 
evidence for the existence of the e-mechanism. 
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